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INTRODUCTION 

The outline shown below describes the material that will be 
covered in the Reliability Testing and Demonstration Lecture IHB: 

RIB. Application of Statistical Methods 

1. Testing with normal units 

2. Determination of confidence limits 

3. Testing with lognormal units 

4. Determination of confidence limits 

5. Testing with bionomial and Poisson evehts 

6. Determination of confidence limits 

A great deal of work has been done by various researchers to 
develop mathematical concepts suitable for reliability studies. The 
interested reader should consult References 1 through 4 for additional 
details pertaining to statistical methods for discrete and continuous 
random variables. 

In these notes effort will be concentrated on four functions: 

(1) Failure, f(t)‘ (2) Reliability, R(t); Failure rate, X; and (4) Hazard 
rate, X’. Since it is usually important to know how well a point esti- 
mate has been defined, some consideration will be given to calcula- 
tion of confidence limits for normal, lognormal and bionomial functions. 
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These notes consider specific cases to show how statistical methods 
can be used in analyzing test data. 

1. Testing with normal units 

A mechanical part is being used where friction, mechanical 
loading and temperature are the principal failure causing stresses. 
Assume that tests to failure have been conducted for these mechanical 
parts resulting in the data shown in table I. 

(a) Calculate the mean-time-between-failures and standard 
deviation. 

(b) What are the hazard rate at 85. 3 K hours and failure rate dur- 
ing the next 10 o 3 K hour interval? 

(c) What are the failure and reliability time functions? 

The mean-time-between-failures and standard deviation can be calcu- 
lated for the data given in table I as follows: 


n 



n 


where 

t mean-time-between-failures 
tf time-to-failure 
n number of observations 
therefore, using the data from table I 

t = I52-K = 75 k hours 


10 


3 


and 



where a = unbiased standard deviation 
2 


n 


t f ] = 7. 50 x 1(T Khr 


>1 


/ n \ 2 2 2 ( 

fy t f ] = (7. 5 x 1(F) = 5.625 x 10* 

\M / 


Therefore, 

„ /57213 - 56250\ 1,/2 A963\ 1/<2 

a= ( 9 ) = (t) = 

(b) The hazard rate, V and failure rate 

follows: 

X 1 - Normal ordinate at 85. 3 K hr 
Normal area 85. 3 K hr to « 

Let Yj = Normal ordinate at 85. 3 K hr 

Zj = Standardized normal variable 
_ t - t _ (85. 3 - 75.0) Khr _ fl Q 
( 7 10. 3 Khr 


(K hr) 2 
10. 3 K hr 

X are calculated as 
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Table 4 (p. 352 of ref. 5) for Z = + 1.0 gives Yj ~ 0.242. The scale 
constant K g for this problem is: 



where 9 = class interval 


Y, = f(t.) = KY, « 1 ° * 1 F X 0. 242 = 2. 35 X NT 4 F/hr 
1 1 s 1 10.3 K hr 


Let R(tj) = Normal area 85. 3 hr to ». From table 3 (p. 351 of 
ref. 5) for Z ^ = +1.0 

Q(tj) - 0. 841 area from -<» to Z j 
Since Q(tj) + R(tj) - 1.000 

R^j) = 1.000 - 0.841= 0.159 


and 


x , _ 2. 35 x 10~ 4 F/hr 
1.59 x 10" 1 




R(t 2 )~ 


= 1 47 x 10'“ Failures/hr 


h - 10. 3 K hr 


Rftg) = Normal area 95.6 K hr to 

z, = (9 5. 6 K- 75.0) K hr = +2 ,o 

* 10. 3 K hr 


From table 3, Q(t 2 ) = 0.977 and R(t 2 ) - 0.023 

\ = ( 1 - 56 <_ 10 4 = 8 31 x 10 -5 Failures / hr 

10. 3 K hr \ 0.159/ 10 3xl0 4 


(c) The constants for f(t) and R(t) are calculated as follows: 


1 - = = 3. 87 x 10 -5 

1.03 X 10 4 x 2.52 

2cj 2 s 2 x (1. 03 x 10 4 ) 2 = 2 . 12 x 10 8 
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Therefore, 

f(t) = 3.87xl0- 5 e-< t - 7 5><104 ) /212X10 8 
R(t) = 3. 87X10' 5 f e"^ ' 7 ‘ 5xl ° 4 ) /2* 12x10® dt 

2. Determination of parameters and confidence limits 

Twenty-five (25) mechanical parts have been tested to failure. 
The mean -time -between failures has been calculated to be 75 K hours 
with cr = 10. 3 K hours (see problem 1). (a) What are the upper and 
lower confidence limits at a 90 percent confidence level ? 

The upper and low^r confidence limits are given by: 

U = F+ V2-T 

Vn 


L = t - K 


where 


t 

K 

a 

n 

Of 

2 


a/2 


a/2 




mean-time-between failures 
normal coefficient 
unbiased standard deviation 
number of samples 

area under one tail 


For the same problem: 
1 - a - 0. 90 


a = 0.10 


- ~ 0.05 


K 


Of/2 


= 1.64 from table 3 of reference 5 


u = 75 K + X — 64 - X 10 - 3 - ~ = 78. 4 K hr 


725 
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L = 75 K - L 04 X 10 ' 3 -ff = 71. 6 K hrs 

725 


This means that 00 percent of the time the mean -time -between - 
failures estimate t for problem 1 with a larger sample size will be 
between 71 600 and 78 400 hours. It should be noted that the sample 
size for problem 1 was only 10 parts. If possible, try to keep n 25 
for estimating normal parameters with the above equations. 

If the sample size, n < 25 then use should be made of the stu- 
dent’s t distribution (see ref. 6). If problem 2 is reworked for a 
smaller sample size of 10, it will be interesting to see the effect that 
sample size has on the size of confidence intervals. 

U = F + V2 ~ 

Vn 

L = F - t a / 2 — 

Jn 


where 

t a / 2 student t coefficient 
s standard deviation 

For this case with a/2 = 0. 05, i/ = n- l = 10-l = 9 

t a / 2 = 2.26 from table IV of reference 6, page 243 

10 


s = 


'57213 


56250 




= 9. 82 K hr 


U = 75 K hr + 


L = 75 K hr - 


2.26 x 9. 82 K hr 


/To 

2.26 x 9.82 K hr 

/ 10 


= 82.0 K hr 


= 68. 0 K hr 
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It will be noted that the smaller sample size gives a larger interval 
of uncertainty for t. 

3. Testing with lognormal units 

A cable used as guy supports for sail experiments in wind tunnel 
testing exhibited the time -to -failure performance data given in table n. 

(a) Write the failure and reliability functions. 

(b) What is the hazard rate at 5710 hoi is? 

(c) What is the failure rate during the next 3000 hours ? 

a. The essential steps for solving this problem are given below: 

30 Table 2 gives the median rank for each ordered position. 

(2) Plot on lognormal probability graph paper (probability x 
2 log cycles) median ranks against failure age as shown in figure 1. 

(3) If a straight line can be fit to these plotted points, then 
the time -to -failure function is lognormal. 

(4) The mean-time -between -failures is calculated by 

t* = log e (t) where t = 6970 hours as shown in figure 1 for a median 
ranks of 50 percent, hence t’ = 8. 84. 

(5) The standard deviation is calculated by 



where t^j = 49 500 hours and t L = 1020 hours 


as shown in figure 1 for a median and [1-rank] of 93.3 percent; hence 



2. 82. Using these constants the expressions 


for f(t) and R(t) can be obtained. 

f(t) = 1- 47x10 -1 e -(t'-8. 84) 2 /l. 59X10 
t 
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The lognormal area from t’ to infinity can be obtained directly 
from figure 1 using the [1-rank] scale. Enter the time -to -failure 
ttf ordinate at 5715 hours; project over to the lognormal life function 
Q(t) and up to the [1-rank] abscissa value of 0.638. Therefore the 
the hazard rate X’ at 5715 hours is: 

X’ = 2. 47x10 — = 3 87x10"^ Failures/hr 
6. 38x10“* 


c. The failure rate during the next 3000 hours is calculated knowing 
the E(tj) = 0. 638 at ttf =5715 hours and obtaining Rftg) = 0. 437 
from figure 1 at ttf = 8715 hours. Therefore, 

X = — - — (l - 5^432. | = l. 05x10”^ Failures/hr 
3xl0 3 \ °* 638 / 
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4. Determination of confidence limits 

It has been shown that the guy supports of problem 3 exhibited a 
reliability of 0. 628 at a ttf of 5715 hours. Consider now the procedure 
for determining the confidence band on this lognormal estimate. The 
d at;*. needed for the graphical construction of the 90 percent confidence 
lines on the lognormal graph of figure 1 is also given in table n. The 
steps necessary to graphically construct the confidence lines in fig- 
ure 1 are as follows: 

(1) Enter the ranks axis with the first 5 percent rank 
value hitting Q(t) the lognormal life function shown in figure 1; 
ordered sample number 3, 5 percei t rank 8. 7. 

(2) Draw a vertical line to intersect Q(t) at point 1 as shown 
in figure 1 . 

(3) Draw a horizontal line to cross the corresponding median 
rank; ordered sample number 3, median rank 25. 9. 

(4) The intersection point (point 2 in fig. 1) of step 3 and the 
median rank line is one point on the 95 percent confidence line. 

(5) Repeat steps 1 through 4 until the desired time -to -failure 
is covered, 5715 hours in this case. 

(6) The 5 percent confidence line is obtained in a similar 
manner. Enter the ranks axis with the 95 percent failure 
rank, 25.9 for ordered sample number i. 

(7) Draw a vertical line which intersects Q(t) at point 3. 

(8) Draw a horizontal line to cross the corresponding median 
rank; ordered sample number 1, median rank 6. 7. 
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(9) The intersection point (point 4 in figure 1) of these two 
lines is one point on the 5 percent confictence line. 

(10) Repeat steps 6 through 9 until the desired time -to -failure 
is covered. 

At 5715 hours the 90 percent confidence interval for Q(t) is from 
figure 1: 19. 7 percent, 69.4 percent. Hence, a 90 percent confidence 
interval for R(t) at 5715 hours is 0. 803 to 0. 306. Incidentally, this 
graphical procedure for finding confidence intervals is completely 
general and can be used on other types of life test diagrams. 

5. Testing with the binomial and Poisson events 

The binomial and Poisson distributions are discrete functions of 
the number of failures which occur rather than time t. A summary 
of these frequency functions is given in figure 1 of reference 4. 

A suspicious lot of explosive bolts is estimated to be 15 percent 
defective due to improper loading density observed in neutron radiog- 
raphy. 

(a) Calculate the probability of one defective unit appearing in a 
flight quantity of four. 

(b) Plot the resulting histogram. 

(c) What is the reliability from the first defect? 

Not much failure density data is available, however, past experience 
with pyrotechnic devices has shown that the binomial distribution applies. 
From the given data: 

q p?r unit number of effectives = 0. 85 
p per unit number of defectives = 0. 15 
n sample size = 4 

Nj possible number of failures = 0, 1, 2, 3, 4 
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The frequency functions corresponding to these constants are given below: 



nl 

(n - N f )! N f ! 


N f n-N f 
P q 


f(N f ) = 


4! 

(4 - N f )! N f ! 


N f 4 -Nr 

p q 


R(N f ) =2 
i=N f 


n! 


(n - j)? V 





pV J 


One easy method to obtain the binomial expansion coefficients is to 
make use of Pascal’s triangle. Pascal found that there was symmetry to 
the coefficient development and explained it as shown in table HI. Col- 
umn 1 gives the sample size n. Column 2 gives the possible number of 
failures. Column 3 gives the binomial expansion coefficients. 

The numbers in the dashed triangle in column 3 are obtained by adding 
the two numbers above the number to get that number; that is, refer to 
dashed insertion the triangle 3+3 = 6. In expanded form f (Nf) becomes 
f(N f ) = q 4 + 4q 3 p + 6q 2 p 2 + 4qp 3 + p 4 
The probability of one defective unit appearing in the flight quantity 
of 4 is given by the second term in the expansion; hence, 

4q 3 p = 4(0. 85) 3 (0. 15) = 0. 37 

The resulting histogram for this distribution is shown in figure 2. The 
probability that 2, 3, or 4 defects will occur as the reliability from the 
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first defect is the sum of the remaining terms in the binomial expansion. 
This probability can be calculated using the following equation: 


R(N f ) = Y i! p j q n 'j 

p2 < 4 - 5) •’ 

However, it is simpler to use the histogram graph and sum the proba- 
bility defects over from 2 to 4. Hence, 

R(2) = 0. 096 + 0. Oil + 0. 001 = 0. 108 
These explosive bolts in their present form are not suitable for use on 
a flight spacecraft as the probability of zero defects is only 0. 522 much 
below the usually desired 0. 999 for pyrotechnic spacecraft devices. 

The Poisson distribution is used to determine the probabilities 
associated with a specified number of failures in the continuum of time. 
Complex electrical components have been shown to follow the Poisson 
distribution. 

Ten space power speed controllers were tested during the Sun- 
flower development program. The time -to -failure test data is given 
in table IV. 

(a) Write the Poisson failure density and reliability functions. 

(b) What is the probability of five failures in 10 000 hours? 

(c) What is the probability that 6, 7, 8, 9, or 10 failures will 
occur or the reliability from the failure ? 

a. Using the data given in table IV , this problem can be solved as follows: 




t=iii 


= -A. 0 86x10 = 8. 59X10 3 hr/failures 
10 
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Hence the Poisson failure density function is: 


3 Nf 


f(N.) = (t/ 8 - 59)< 10 ) e -t/8. 59x10 


N ? 


and the reliability function is: 
10 


R(N f ) <L/ 8 ' SSxlO 3 ) 3 a -t/8. 59X10 3 


j.' 


3=1 


b. The probability of five failures f(5) in 10 000 hours makes use of 

the ratio (t/F). Using this ratio, f(5) becomes 

t_ _ l.QxlQ 4 _ 1 16 
* 8. 56x10 3 

f(5) = ill 16 ^ 5 e = 2 -.09 x Q. 314 = 5 47x10 -3 
5! 1.2X10 2 

One easy method to calculate the term (1. 16) 5 is as follows: 
log (1. 16) 5 = 5 log 1. 16 = 5(0. 148) = 0. 740 
(1.16) 5 = 2.09 

c. The reliability from the 5 th to the 10 th failure is the sum of the 

remaining terms in the Poisson expansion. This probability can 
be calculated using the following equation 


10 

R(N f )=£ 

3=6 


0.314 (1. 16) j 
j! 


R(6) = 0.0013 
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6. Determination of confidence limits 

When an estimate is made using discrete distributions, it is ex- 
pected that additional estimates of the same parameter will be close to 
the original estimate. It is desirable to be able to determine upper and 
lower confidence limits at some stated confidence level for discrete dis- 
tribution estimates. The analytical procedure for determining these 
intervals is simplified by using specially prepared tables and graphs. 
Useful tables for the binomial distribution are given in references 5, 8, 

9, and 10. 

A prior calculation showed that the probability of one defective pyro- 
technic unit appearing in a flight quantity of four was 0. 37. 

What are the upper and lower confidence limits on this estimate at 
a 90 percent confidence level? 

If the number of defectives is r and the confidence level is y, this 
problem has the constraints listed below: 

n = 4 r = 1 y = 90 percent 

Using these constraints, the upper U and lower L confidence limits 
can be obtained from table I in reference 8. 

L = 0.026 
U= 0.680 

This means that with a 90 percent confidence the probability of one 
defective bolt appearing in a flight quantity of four is in the interval from 
0.026 to 0.680. 

The reliability from the 5th to the 10th failure for speed controllers 
was found to be 0. 0013 in a previous problem. What are the upper and 
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lower confidence limits on this estimate at a 95 percent confidence level? 

The variation in t can be found by using Chart I, page 23 from ref- 
erence 10. Enter Chart I on the 5 percent line at the left hand end of the 
5 interval, here T/t^ = 10. 5; then tj = 10 t/T/tj = 8. 57x10^/10. 5 = 
,8160 hours. Using the left hand end of the 4 interval T/fg = 9. 25; then 
tg = 8. 57xi0 4 /9. 25 = 9530 hours. One easy method to find Q(6) is to use 
figure 6-1 of reference 5, page 61. The t/t ratios of interest are 1.05, 
1. 16, and 1.22, respectively. For these ratios with Nj = 5, the values 
of Q(6) from figure 6-1 are 0. 997, 0. 9987, and 0. 9992, respectively. 
Since the sum of the last five terms is desired, R(6) is 0. 003, 0. 0013, 
and 0. 0008, respectively. 

This means that the probability of the 5th to the 10th failure of a 
speed control occurring is in the interval from 0. 0008 to 0. 003 at a 
confidence level of 95 percent. 



16 

REFERENCES 

1. Anon. : Reliability Theory and Practice. ARINC Res. Corp. , 

Washington, D.C., 1962. 

2. Anon. : Reliability by Design. General Electric Co. , Defense 

Elec. Div. , Waynesboro, Va. , 1964. 

3. Bazovsky, Igor: Reliability Theory and Practice. Prentice -Hall, 

Inc. , 1961 

4. Lalli, V. R. : "Reliability Testing and Demonstration, " presented 

at the Sixth Annual Reliability Engineering and Management 
Institute, Univ. Arizona, Tucson, Ariz., Nov. 5-7, 1968. 

5. Calabro, S.R. : Reliability Principles and Practices. McGraw-Hill 

Book Co. , Inc. , 1962. 

6. Hoel, Paul G. : Elementary Statistics. John Wiley & Sons, Inc. , 

1960. 

7. Anon : Failure Distribution Analyses Study. Vols. I, n, and HI. 

Computer Applications, Inc. , Aug. , 1964. 

8. Lochner, R.H. : Estimation and Prediction Using the Binomial Dis- 

tribution. Reliability Res. and Ed. Dept. , General Motors Corp. , 
Milwaukee, Wise. , 1963. 

9. Leone, F. C. , et al. : "Percentiles of the Binomial Distribution, " 

Case Inst. Tech. , Cleveland, Ohio, 1967. 

10. Lochner, R. H. : Reliability Calculations for Exponential Population. 
Reliability Res. and Ed. Dept. , General Motors Corp. , Milwaukee, 
Wise. , 1963. 



1? 


TABLE I. - TEST DATA FOR A MECHANICAL PART 


Ordered sample number 

f f , 

r 2 

''’f > 


K hr 

(K hr) 2 

1 

60 

3600 

2 

65 

4225 

3 

68 

4624 

4 

70 

4900 

5 

75 

5625 

6 

75 

5625 

7 

80 

6400 

8 

83 

6889 

9 

85 

7225 

n = 10 

90 

8100 

Totals 

750 

57 213 


TABLE H. - CABLE TIME -TO -FAILURE DATA 


Ordered sample 
number 

Time -to - 
failure, 
hr 

Median 

rank 3. 

5% Rank 1 

95% Rank 3 - 

1 

1 10C 

6,7 

0. 5 

25.9 

2 

1 890 

16. 2 

3.7 

39.4 

3 

2 920 

25.9 

8.7 

50. 7 

4 

4 100 

35.5 

15.0 

60.7 

5 

5 715 

45 2 

22.2 

69.7 

6 

8 720 

54.8 

30.3 

77.8 

7 

12 000 

64.5 

39.3 

85.0 

8 

17 500 

74.1 

49.3 

91.3 

9 

23 900 

83.3 

60. 6 

96.3 

n a 10 

46 020 

93.3 

74.1 

99.5 


a From tables 2, 5, and 15 o i reference 7 . 
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TABLE HI. - PASCALS TRIANGLE 


FOR BINOMIAL COEFFICIENTS 


Sample 

Possible 

Binomial 

size 

failure 

coefficients 

1 

2 

1 

2 

3 

1 2 1 

3 

4 

1 \3 3/ 1 

S3 

II 

5 

1 4 \6/ 4 1 


TABLE IV. - SPEED CONTROLLER 


TIME -TO -FAILURE DATA 


Ordered sample number 

T im e -to -f ailu re , 
hr 

1 

3 520.0 

2 

4 67L2 

3 

6 729, 3 

4 

7 010.0 

5 

8 510.2 

6 

9 250. 1 

7 

10 910.0 

8 

11 220. 5 

9 

11 815.6 

10 

12 226.4 

Total 

85 866.3 
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Figure 1. - Lognormal life test diagram. 



NUMBER OF FAILURES, N f 
Figure 2, - Explosive bolts histogram. 




